Derivation of the Gauge Link in Light Cone Gauge 

Jian-Hua Gaqj 

Department of Modern Physics, University of Science and Technology of China, 

Hefei, Anhui 230026, China and 
Q ■ Department of Physics, Shandong University, Jinan, Shandong 250100, China 

(n: 

j^. (Dated: September 1, 2010) 

<^ ■ Abstract 



m 



In light cone gauge, a gauge link at light cone infinity is necessary for transverse momentum- 
dependent parton distribution to restore the gauge invariance in some specific boundary conditions. 
Qh' We derive such transverse gauge link in a more regular and general method. We find the gauge link 






X 



at light cone infinity naturally arises from the contribution of the pinched poles: one is from the 
quark propagator and the other is hidden in the gauge vector field in light cone gauge. Actually, 



^^ \ in the amplitude level, we have obtained a more general gauge link over the hypersurface at light 



cone infinity which is beyond the transverse direction. The difference of such gauge link between 



in 

o 

.^ . semi-inclusive deep inelastic scattering and Drell-Yan processes can also be obtained directly and 

•Q \ clearly in our derivation. 
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I. INTRODUCTION 

Nucleon structure functions are physical observables and can be measured in deep inelastic 
scattering (DIS). In the naive parton modelfll], the structure functions are expressed in terms 
of the probability of finding quarks and gluons in the parent nucleon. In collinear QCD 
factorization formulas, such structure functions can be given by compact operator matrix 
elements of the target [2] 

q{x) = ^ /"^e-«P^^"(P|^(y-,Ox)^£b-,0^;0,Oj^(0,0^)|P), (1) 



where 



C[y',6r,0,6^] = Pexpl-tgJ' dCA+iCO^)) , (2) 



is the gauge link between the quark fields, which arises from final state interactions between 
the struck quark and the target spectators. In Eq. ([1]) and Eq. ([2]), all fields are evaluated at 
equal y~^ = 0. Since structure functions, as physical observables, should not be dependent 
on the gauge that we choose, it is necessary to introduce such gauge link to ensure the 
gauge invariance of matrix element. In the light cone gauge A"*" = 0, where the path-ordered 
exponential in Eq. ([2]) reduces to unity, we can identify the quark distribution in Eq. ([T]) as 
a probability distribution as we made in naive parton model. Actually, in collinear structure 
function such as Eq. ([1]), we can always select a clever gauge to vanish the gauge link. 
But when we consider the transverse-momentum dependent quark distribution, such naive 
manipulation will result in inconsistency. In the nonsingular gauge, in which the gauge 
potential vanishes at the space-time infinity, the transverse-momentum parton distribution 



is defined in the literature as 



y 



g(^,^±) = l j '^-f^e-^^y-^^^^y^ 



2 J 271 (27r)2 

X (P|V^(y-, yA.)iO[oo, y^- y-, f J£[oo, 0^; 0, J^(0, 0^) |P) , (3) 

where 

C[oo,yr,y',y±\ = Pexpi-ig dC A-^ (C , y±) j , (4) 

and all fields are evaluated at equal y^ = 0. From Lorentz invariance, parity invariance and 
time reversal invariance, the transverse- momentum parton distribution can be decomposed 



into the following expressions, 

q{x,k^) = f{x,k^) + S-{fxk^)f^^{x,k^)/M (5) 

where S is the spin of the target nucleon and p is a unit vector along the direction of the 
target momentum in infinite momentum frame. The function f^{x, k±) is just the Sivers 
function and can contribute to single spin asymmetries. It is verified in Ref. 5] that the 
Sivers function vanishes unless there is the gauge link in Eq.(jll), which is yielded by the 



final state interactions 



0- 



In the light cone gauge, however, it seems as if the gauge link 
in Eq. (j3]) would become unity and the final interaction vanish accordingly too. Hence 
there will be inconsistent results from different gauges, which is impossible, since physical 
observables should not depend on the gauge by choice. Ji and Yuan in [7[ have shown that 
the final state interaction effects in single spin asymmetry can be recovered properly in the 
ight cone gauge by taking into account a transverse gauge link at y" = +oo. Further in 
Belitsky, Ji, and Yuan demonstrate the existence of extra leading twist contributions 
from transverse components of the gauge potential at the light cone infinity. It turns out 
that these contributions just form a transverse gauge link in light cone gauge. In this paper, 
we will give another more regular and systematic method to obtain such transverse gauge 
link in light cone gauge. We find the gauge link at light cone infinity will arise naturally 
from the pinched poles, one of which is provided by the quark propagator and the other is 
hidden in the gauge vector field in light cone gauge. Actually, it turn out that we obtain a 
more general gauge link over hypersurface y~ = oo, instead of only transverse gauge link. 
The difference of such gauge link between semi-inclusive deep inelastic scattering (SIDIS) 
and Drell-Yan (DY) processes can also be shown directly and clearly in our derivation. The 
paper is organized as follows: in the next section, we will introduce some kinetics definitions 
and notations which will be involved all through our paper. In Sec. Ill, we would like to give 
a brief review on the singularity in light cone gauge and different prescriptions for different 
light cone pole structures. Then in Sec. IV, we will devote to deriving the gauge link in light 
cone gauge in SIDIS process. In Sec.V, we will deal with the DY process and compare it 
with the SIDIS process. A very short summary is given in the end. Other relevant work on 
the transverse gauge hnk can be found in the literature [ol, [lO| ■ 



II. SOME DEFINITIONS AND NOTATIONS 

In studying SIDIS or DY process, it is convenient to choose the hght cone coordinate 
system in which we introduce two hghthke vectors n^ and n'^, 

n^ = (0,1,0^), n'^ = (1, 0, 0^) , n-n = l. (6) 

With these basis vectors, we may write any vector k^ as (A;^, A;~, k±), where k~^ = k-n, k~ = 
k-n. For example, in SIDIS process, we choose the proton infinite momentum frame, in which 
the proton's momentum and the virtual photon's momentum are given by, respectively, 

pt" = p+n^^ qf" = -xbP^ + — ^ — n'^. (7) 

2xbP^ 

where Xb = Q"^ /'^V " ^ and Q^ = —q^. 

In order to make the derivation more compact and elegant in the following sections, let 
us introduce some notations. For any momentum vector k^ and the gauge potential vector 
A^, we will manipulate the following decomposition: 

F = F + xp^. A" = A^' + A^vf (8) 

~ — # ~ — * 

where k^ = {0,k~ ,kj_), x = k~^/p~^, and A^ = (0, A^ , A±). For any coordinate vector y^ , 
we will make the dual decomposition, 

y^^ = r + y"n^ (9) 

where y'^ = {y~^, 0,y±). When there is no confusion, we will rewrite y'^ as {y~, y). With such 
notations, we have k ■ y = k ■ y + xp^y~ , and in light cone gauge where A'^ = ,we also 
have A^ = A^. It should be noted that in light cone coordinate, the covariant vector and 
contravariant vector are related by A~^ = y4_, A~ = A^ and A-^ = —A±. 

III. SPURIOUS SINGULARITY IN LIGHT CONE GAUGE 



The light cone gauge n • A = is widely used in perturbative QCD calculations [Uj, ll2| , 
and under such a physical gauge condition, the probability interpretation is expected to 
hold. The Yang- Mills theories, quantized in light cone gauge, have been studied by several 



authors 



13 



14| . However, when we calculate with the gauge propagator in such gauge 



in perturbation theory, we have to introduce some spurious pole to regularize associated 
light cone singularity. There have been a variety of prescriptions suggested to handle the 
singularities 15l-ll9| , in which most attempts were pragmatic. The literature [sl, |20| states 



that in general, in light cone gauge, the gauge potential can not be arbitrarily set to vanish 
at the infinity, the spurious singularities, characteristic of all the axial gauges, are physically 
related to the boundary conditions that one can impose on the potentials at the infinity. In 
our paper, we will consider three different boundary conditions as in |8|, i.e. 

Advanced : A{oo,y) = 

Retarded : A{—oo, y) = 

Antisymmetric : A{—oo,y) + A{oo,y) = 0. (10) 

The typical integration we will meet with in our derivation is the Fourier transformation of 
the gauge potential such as, 

/oo 
dy~e''^y-Apiy-,y) (11) 

'OO 

Manipulating this integration by parts, we obtain 

dy-e^'^y-Ap{y-, y) = [—] / dye^'^y' d^A^{y', y) (12) 

-oo "^ J —oo 

where d~^ = d- = d/dy~ . Since the boundary condition is set, the term [pr] can be 
regularized by definite prescription. 





3tric : 


I 

^k+^ 
i , 

i 

^k+^ 


I 


Retarded : 


k+ — ie 
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k+ + ie 




2^k+ + ie 



where the last propose is just the conventional principal value regulation when the antisym- 
metry boundary condition is assigned. Hence, we notice that there is a secret pole structure 
in gauge potential in momentum space. We will show that it is just this pole that will 
contribute to the final gauge link at the light cone infinity. 

The easiest way to illustrate the validity of such regularization is just to set 

Advanced : 9{—y~) 

Ap{y~) = I Retarted : e{y-) (14) 

Antysymmetry : i [0{y') - 9i-y')] 



where the function 9{y~) is the usual step function. It is a trivial exercise to show that they 
can result in the proper pole structure as we present in Eq. f lTOl) . 

As we mentioned above, in the light cone gauge, we can not impose on the gauge potential 
the boundary condition both y4p(+oo, y) = and y4p(— oo, ij) = 0. We can only choose either 
of them as the boundary condition to remove the residual gauge freedom and the other one 
will be subjected to satisfy the field equation or the request that the total gauge energy 
momentum is finite. However, as a matter of fact, we can still impose a weaker condition, 
that the gauge potential must be a pure gauge. In the Abelian case, 

ip(±oo, y) = dp(j){±oo, y) (15) 

or in the non- Abelian case 

ip(±oo, y) = a;"^(±oo, y)dpu{±oo, y) (16) 

where u = exp(i0). In the non- Abelian case, Ap = A'^t"' and = (j)°'t"' where t"" are the 
generators of non-Abelian group in the fundamental representation. Keeping the leading 
term in the Tailor expansion of u around (p, we recover the same expression as Eq. ( ITSll in 
the Abelian case. It follows that 



(j){+oo,y) = - d^-Ai+oo,0 (17) 

where the integral runs over any path on the hypersurface y^ = oo. Notice that this equation 
always holds for Abelian gauge potential, and holds for the non-Abelian case only when the 
is small. It will be interesting thing to investigate what the nonleading terms contribute 
to in the non-Abelian case, which is beyond the scope of this paper. We will show that the 
linear term, such as in Eq. ( 1T7|1 will lead to the gauge link at the light cone infinity. 



IV. GAUGE LINK IN LIGHT CONE GAUGE IN SIDIS 



In DIS process, the hadronic tensor is defined by 



^'^ = ^E/^(2^)''^^'^(^x+P.-p-?)(^|jnO)bj,X)(p,,X|/(0)|P). (18) 




k+q 



^ 



FIG. 1: The tree diagram in DIS process 




FIG. 2: The one-gluon exchange diagram in DIS process 



The tree scattering amplitude corresponding to Fig. [T] reads 



M^ = (pj,X|/(0)|P)(o) = uik + q)r{XmO)\P) , 



(19) 



where k denotes the momentum of intial quark scattered by the photon with momentum q. 
The one-gluon amplitude in light cone gauge corresponding to Fig. |2] reads, 



Mf 



d% 



(2vr)4 
xu{k + q)Y^ 






{ki + qY + ie 
The quark propagator can be decomposed into two parts 



{X\A,,{y^)rHO)\P) 



h 



h+ 



(xi — xi + ie 



+ i) 



(20) 



(21) 



(/ci + qf + ie 2p- Cki + q) 

where ki = (xip"*", k^ , kij_) with xi = k^/p'^ = xb + k'j_/2p ■ {ki + q) is determined by the 
on-shell condition (/ci + g)^ = . Actually, to obtain the Eq. ( 12T|) . we have neglected the 



contribution 

e(-A;r - q- 



2p ■ (fci + q) 



^1+^ ..^ 



(22) 



[Xi — Xi — It) 

which will contribute at higher twist level since they vanish in the limit q~ — )■ +00. The 
last term in Eq. ( pT]) is the so-called "contact" term of a normal propagator which does not 



propagate along the light cone coordinate 2l|. Such a contact term will always result in 



higher twist contribution and does not contribute to gauge link at all. Hence, when we are 
considering the leading twist contribution in our following derivation, we can just drop such 
contact terms and only keep the pole terms, i.e. the first term in Eq. ( 1211) : 



^^^1 f „,3 • [ P'^dXi I i(^k-h)-yi+i{x-x,)p+y- 



where another notation Mf with an extra " is introduced to remind us that the only pole 
term is kept, and we have also separate the integral over xi and yf from the others which 
means we will finish integrating them out first in the following. Before proceeding further, we 
should first choose a specific boundary condition for the gauge potential Ap at infinity. Let 
us start with the retarded boundary condition A{oo,y) = 0. Using the Eq. (iT2l) accordingly 
which corresponds to retarded boundary condition, we have 

^^1 f ^3- f dXi^ f r- i(k-ki)-yi+i(x-Xi)p+y- 






2p ■ {ki + q) {xi -Xi+ it) (x - Xi + ie) 
x(X|9+i,,(yi)7^^(0)|P). (24) 



Now we can finish integrating over Xi and Hi first, 

'^'^1 / 7 - i{x-x-i)p+y- ^ * 



dy^e'^^-^^^^-y :-—-- d^ Ap,{y,) 

ZTT J (xi — xi + le) [x — Xi + le) 



1 



1 

X — Xi 

1 

X — Xi 



dy^ (e{y')e'^^-'^^^^y- +e{-y-)) 9+ip,(yi) 
dyi {e{y-) + 0{-y')) d+Ap,{yi) + higher twist 
Ap-^ (+C)0, yi) + higher twist , (25) 



where only the leading term in the Tailor expansion of the phase factor e*(^~^i)P ^ is kept, 
because the other terms are proportional to (x — Xi)" = [k']_/2p ■ {k + q) — klj_/2p ■ {ki + g)]" 
(n > 1), which will contribute at higher twist level. Only keep leading twist contribution 
and inserting Eq. ( 125|) into Eq. i HM . we have 



xu(k + q)Y' *^t* .. '- (A1/i,.(+«3.i/iW>(0)|P). (26) 

2p ■ [ki + q)x — Xi 

Using Eq. ( !T5|) and performing the integration by parts over iji where dp — )■ —i{k — ki)p, we 
obtain 



^^ = /^/^'^^ ^''"''^ 



x^(A; + q)C^- |^,)_fcti_^l_(x|0(+oo, mmO)\P) . (27) 

2p ■ (fci + g) a; - xi 

To carry out the matrix algebra further, we note that 

~^-h = {^ + ^)-{h + ^)- (x-x,)^, (28) 

together with the on-shell conditions 

u{k + q){^ + ^) = 0, and {!:, + ^)^ = . (29) 

Using these equations, we reduce the Mi into 

yiu{k + q)m+ i)—^ -(^|0(+oo, yi)V^(O)lP) 

2p ■ (/ci + g) 

= M(A; + g)(X|z0(+oo,O)V'(O)|P) 

xuik + q)§- m ^ -(^l0(+oo, 2/i)^(0)|P) . (30) 

2p ■ (/ci + g) 

Since the last term in Eq. (!30|) only contribute to higher twist, keeping only the leading 

twist contribution, we finally obtain. 

Ml = u{k + q){X\i(P{+oo,0)ilj{0)\P) . (31) 



So far, the previous derivations have been restricted to the retarded boundary condition 
where A{—oo, y) = 0, now let us turn to the other two boundary conditions. When we 
assign the advanced boundary condition A{+oo, y) = 0, which means that we should choose 
the advanced one in Eq. flT^ and Eq. flT^ . Such a sign change in the pole structure will 
lead to replacing the integration in Eq. (12^ by, 



p [dy^e^^-^^^^'y-- ^^—— '—-d^A,,{y,) 

ZTT J [Xi — Xi + le) [X — Xi — le) 

= -^ [dy^ fe{y-)e^^^-'^^P^y~ - Oiy-)) 9+i,,(yi) 

It/ 

= higher twist . (32) 

We note that, different from retarded case, the leading contributions from two poles have 
canceled each other completely, and there will be no gauge link at all. As shown by [8], all 
final state interactions have been included into the initial state light cone wave functions. 

If we choose the antisymmetry boundary condition, which corresponds to the principal 
value regularization, we have 

d^l I _7_ - J{x-xi)p+y- l 



-^ / dy^e^^^-^^^^^y-- ^— -rPVy ,d^A,,{y,) 

2tt J [xi - Xi + le) [x - Xi) 

^ ^dy^^ {2e{y-)e'^^-'^^^^y- - e{y') + 0{-y-)) d+A,,{y^) 



X — Xi J 2 

1 /■ 1 
^— - / dyi- {e{y-) + e{-y-)) d+Ap,{yi) + higher twist 

1 



Ap (+00, yi) + higher twist , (33) 

X — Xi 

where PV denotes principal value. The above result appear the same as the one in the 
retarded boundary condition. The difference between retarded and principal value regular- 
ization is that final state scattering effects appear only through the gauge link in principal 
regularization, while they appear through both the gauge link and initial light cone wave 
functions in retarded regularization. Such detailed discussion and illustration can be found 
in Ref.|8|. In the above derivation, we notice that the pinched poles are needed to pick up 
the gauge potential at the light cone infinity, which will be shown to result in the gauge 
link that we expect. In the following, we will only concentrate on the retarded boundary 
condition in the following derivation. 




FIG. 3: The two-gluon exchange diagram in DIS process 



Now let us consider further the two-gluon exchange scattering amplitude in Fig. El 

-{X\A,,{y2)A,,{y{)ri,mP) . (34) 



M^ 



d'^k2 d'^ki 
(27r)4 (27r)4 

xu{k + q)Y \, ^.^\t . . Y' — 



{k2 + qY + ie {ki + q^ + ie 
Just following what we did in the Mf, we drop the contact terms which do not contribute 
in leading twist level and label the residual terms as M2, which is given by 



M^ 



d^k2 d^ki f ,0 f p~^dx2P~^dxi 
d y2d yi 



(27r)3 (27r)3 J " """ "" J 27r 2tt 

y^^i{k-k2)-y2+i{k2-ki)-yi+i{x-X2)p+y2+i{x2-xi)p+yl 



dy^ dy^ 



P2_ 



%+il 



-7 



Pi. 



k+^ 



xu{k + g)7 

2p ■ {k2 + q) ' 2p- {ki + q) [^2 -X2 + it) [xi -Xi + le) 

x(X|i,,(i/2)ip,(t/i)7^^(0)|P). (35) 



Still with the help of the regularization in Eqs. (1121) and ( 11311 . let us do integrating over x\ 
and yf first, 



M^ 



,3. ,3- / dx2P^dxi , 
d y2d yi I -77- ^ / dy2 dy^ 



d^k2 d^ki 

(27r)3 (27r)3 J " """ "' J 27r 27r 



P2 



^2+^ 



-7 



Pi. 



1^1- 



xu(k + g)7 . , . , ..,/-- • \ / -- • N 

2p ■ {k2 + q) 2p- {ki + q) {x2 - X2 + le) {x-X2 + ^e) (xi - Xi + le) 

x{X\d+A,,{y2)Ap,iy^)YH0)\P) 

d k2 d ki [ jZ„;_Ai„; fP dXi^ f ,^ ^i^k^k2).y2+i{k2-k^yn+i{x-xi)p+y- 



(27r)3 (27r)3 



P2 



d%2d%i 

k+i 



-7 



Pi. 



27r 



x-u(A; + g)7 . , . / ^ • n / - • n 

2p ■ (^2 + g) 2p ■ (fci + g) (a; - 0:2 + ^e) (xi - xi + le) 

x(X|4,(+oo,y2)ip,(l/i)7"^(0)|P). 



(36) 



Further integrating over x^ and y^ , which is totally the same as what we did with xi and 
2/j~. The results read 



2p ■ {k2 + q) 2p- (fci + q) (x - X2 + ie) {x - Xi + ie) 

X {x\A,, (+00, y2)ipi (+00, m)Y^mp) 

d h d ki f 3 3 j(fc_fc2).y2+i{fc2-fci)-yi 



(27r)3 (27r)3 



d-'md-'yi e' 



Xu{k + q)Y^ fe+^ 7" ^^+^ 



2p ■ (A;2 + g) 2p ■ (fci + q)ix-X2 + ie) {x - Xi + ie) 
x{X\d,,^{+oo,y2)dp,(P{+oo,m)rH0)\P) ■ (37) 

Now we are in a position to perform integrating over /c2 and 1/2. Thanks to the integration 
by parts and the algebras given in Eq. (!28|) and Eq. ( 129|) . we obtain 



MM ^ / f^ f ^Sy^ e<'-~'^>y^u{k + q)Y' ^^^ ^ ^ 



(27r)3 y 2p ■ (fci + q)ix-xi + ie) 

X (XK0(+oo, yi)5p,0(+oo, yi)Y^{0)\P) 



rf=^2/i e^('=2-'^^)-*^M(A; + g)7' 



(27r)3 J 2p ■ {ki + q){x-xi + 



le) 



x(X|^9p,02(+oo,yi)7'^V^(O)|P) (38) 

Repeat what we did with /c2 and 1/2 above, and we can finish integrating over ki and y\ and 
finally arrive 

M^ = uik + q){X\'-<p'i+oo,0)rHO)\P) (39) 

All through calculating M2, as we did with Mi, we have neglected the higher twist con- 
tributions and only keep the leading twist terms. From Mi to M2, it is obvious that our 
procedure can be easily extended to n-gluon exchange amplitude M„ in Fig. HI which is 
given by 




^ 



k,+q 


k^+q 
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FIG. 4: The n-gluon exchange diagram in DIS process 



M„ 



rf^A;, 



n___i_/73 • „«{A;„-fc„-i)-i/n+J(fcn-i-fcn-2)-yn-i+ +i{k2-ki)-y2 

TT i^ i(^y-g«(a;n+l-3::n)p+?/n+i(3::n-a;„-i)p+y^_i+ +iix2-3:i)p+ y^ 



i=i 



27r 



x-u(A; + q)y 



Pn I n 



k+i 



7 



pi 



h+^ 



2p ■ {K + q) 2p • (/ci + q) 

1 1 



(x„ — Xn + ie) [Xi — xi + ie) 



(40) 



We first finish integrating from XniUn ^o Xi,yi one by one. Keeping the leading twist 
contribution, we have, 



M„ 



J^ H^k- - - - - 

T-j- _^^_|\7_ 73 • ^i(fc„+i-fc„)-y„+i(fc„-fc„_i)-y„_i+ +i{k2~ki)-y2 

11 (2^)3^ ^^-^ 

7 = 1 ^ ' 



xu{k + q)Y" 



k+i 



2p ■ {kn + q) 



7' 



pi. 



k+i 



2p ■ (/ci + q) 



X- 



l^^n+l ~ Xn) [X2 — a^ij 

x{X\A,^i+oo,yn)Apn-ii+oo,yn-i) A,,{+oo,yi)ij{0)\P) , (41) 



or using 



Apni + ^^Vn) = (9p„0( + OO,y„), 



(42) 



we rewrite it as 



JWn - y 11 (27r)3 ^^ 



X 



2p ■ (/c„ + g) 2p ■ (fci + g) 

1 



(£„+i-x„) (X2-Xi) 

x(X|ap„0(+oo,yO9p„_,0(+oo,y„_i) ap,0(+oo,yi)^(O)|P). (43) 

Continue to integrating over from kn and ijn to k\ and ?/i one by one, we can finally have 

M„ = n(A; + g)(X|^0"(+oo,O)^(O)|P). (44) 

As a final step, we should resum to all orders and obtain 

J2^n = u{k + q){X\exp{i<f){+oo,0))tP{0)\P), (45) 

n=0 

or the more conventional form 

^M„ = M(A; + g)(X|Pexp -2 / d^ ■ i(+cx),e) W'(0)|P) , (46) 



ra=0 



where Pexp f — i J.°° d^ ■ A{+oo,^)] is just the gauge link that we tried to derive. It should 
be noted that the gauge link we obtain in the final result Eq.( l46l) is over the hypersurface 
at light cone infinity along any path integral, not restricted along the transverse direction, 
which means that it is more general than what Belitsky, Ji and Yuan have obtained in 



s 



Ref.|8| 



V. GAUGE LINK IN LIGHT CONE GAUGE IN DY 

Now let us turn to the DY process, which is represented in Fig. (|5]), where, for brevity, 
we have fixed the target to be a nucleon and the projectile to be just an antiquark, q is 
the virtual photon's momentum and q — k and p is momentum of the projectile and target 
respectively. Such simplifying does not lose any generality when we are only considering how 
to derive the gauge link, but it will be more convenient and manifest to compare with the 
SIDIS process. We still choose the light cone coordinate system, and use the two lightlike 




FIG. 5: The one-gluon exchange diagram in DY process 

vectors n^ and n^ to fix "plus" and "minus" directions. All the notations and conventions 
are the same as in DIS. The one-gluon exchange amplitude reads 

xuiq - k)Y' ^^j~J\^^ {X\A,,iy,)ri;mP) . (47) 

Dropping the contact terms and assigning the retarded boundary condition, we rewrite it as 



pi. 



xu{q - k)y . , ^ ■ ^ , . x 

2p ■ {ki - q) [xi - xi - le) [x - Xi + it) 

x(X|9+4,(yi)7'^^(0)|P). (48) 



It should be noticed the difference of the pole structure between Eq. fHHj) and Eq.f l24j) . Just 
like we did in the SIDIS, we can finish integrating over xi and y^ first, 

"^^^ ^ dy^e'^^-^-'^^'-y-- \ ^ d^Ap,{y^) 



27r J (xi — Xi — ie) [x — Xi+ ie) 






= dy^ {0{-y-) - e{-y-)) d+Ap,{y^) + higher twist 

X — Xi J 

= higher twist . (49) 

Opposite to the case in SIDIS, the retarded boundary condition does not lead to the gauge 
link at the light cone infinity and hence all the final state interaction effects must be shifted 



into the initial light cone wave functions. For the advanced boundary condition, we have 

/V^ _ f d^h /• 3 f ^ f ^^ iCk'kiyyii{x-xi)p+y- 



2p ■ {ki - q) {xi -Xi- it) [x-xi- it) 
X (X|9+i,, (1/1)7^^(0) |P). (50) 



Finish integrating over Xi and y^ : 



"^^^ ' rfyre^("-"i)*'^^"- "- — -d^Ap,{yi) 



2tt J [Xi — Xi — it) [x — Xi — it) 






-^— / dy^ {e{-y ) + e{y )) d+ Ap,{yi) + higher twist 

Apj (—00, y\) + higher twist . (51) 



X — X\ 
If we choose the antisymmetry boundary condition, we have. 



/ ^ / dy-.e^^-^^^^^y- ^ -VN-^-d^A,, (y,) 

J ZlT J [Xi — Xi — tt) [X — Xi) 

dy^\ (^-2^(-y-)e^(-^0P^^- - e{y-) + e{-y-)) d+ A,,{y,] 



X — Xi 



X — 

1 



^— I dVil (^(^~) + ^(-^~)) 9'^^pi(y^) + higher twist 
Xi J Z 



-Ap-^{— 00, yi) + higher twist . (52) 

The difference between advanced and principal value regularization in the DY process is 
that final state scattering effects appear only through the gauge link in principal value 
regularization, while they appear through both the gauge link and initial light cone wave 
functions in advanced regularization. It follows that, 

MiiDY) = u{q-k){X\ i-^ f di-A{-^,i)\^mP)- (53) 

Just as in the SIDIS process, the pinched poles are indispensable to produce a finite contri- 
bution in the leading twist level. Following the same line as the one carried out in SIDIS, 
we show that the gauge link in advanced or antisymmetry boundary conditions is given by 

Y,Mn(DY) = u{q-k){X\Pe^Y>{-^ f rf^ • i(-oo,o) V^(0)|P) . (54) 



n=0 



It should be noted that the hght cone infinity y~ = +00 has been replaced by y~ = —00, 
refiecting that the gauge link arises from the initial state interactions rather than from the 
final state. 

To summarize, in light cone gauge, we should choose a specific boundary condition first 
to fix the residual gauge freedom. Using the proper regularization corresponding to specific 
boundary condition, we can obtain the residual gauge link at infinity along the light cone 
coordinate. We find the gauge link at light cone infinity arises naturally from the pinched 
poles: one is from the quark propagator and the other is hidden in the gauge vector field 
in light cone gauge. Actually, it turns out that we obtain a more general gauge link over 
hypersurface y~ = 00, which is beyond the transverse gauge link. The difference of such 
gauge link between SIDIS and DY processes can also be obtained directly and clearly in our 
derivation. We expect our regularization method will also be valuable to make it possible 
to perform higher twist calculations in light cone gauge more unambiguously. 
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